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Production of small-scale Alfven waves by ionospheric depletion, 
nonlinear magnetosphere-ionosphere coupling and phase mixing 

A. J. B. Russell/ A. N. Wright,^ and A. V. Streltsov^ 

Abstract. Rockets and satellites have previously observed small-scale Alfven waves in- 
side large-scale downward field-aligned currents and numerical simulations have associ- 
ated their formation with self-consistent magnetosphere-ionosphere coupling. The origin 
of these waves was previously attributed to ionospheric feedback instability, however we 
show that they arise in numerical experiments in which the instability is excluded. A 
new interpretation is proposed in which strong ionospheric depletion and associated cur- 
rent broadening (a nonlinear steepening/wavebreaking process) form magnetosphere-ionosphere 
waves inside a downward current region and these oscillations drive upgoing inertial Alfven 
waves in the overlying plasma. The resulting waves are governed by characteristic pe- 
riods, which are a good match to previously observed periods for reasonable assumed 
conditions. Meanwhile, wavelengths perpendicular to the magnetic field initially map to 
an ionospheric scale comparable to the electron inertial length for the low-altitude mag- 
netosphere, but become shorter with time due to frequency-based phase mixing of bound- 
ary waves (a new manifestation of phase mixing). Under suitable conditions, these could 
act as seeds for the ionospheric feedback instability. 



1. Introduction 

Small-scale Alfven waves are regularly observed inside 
large-scale magnetospheric current systems [e.g. review by 
Stasiewicz et ai, 2000, and references therein] and several 
studies have shown an association with downward current 
regions [Karlsson et al, 2004; Keiling et ai, 2005; Johansson 
et al, 2004; Wright et ai, 2008, e.g.]. They are speculated to 
play a role in electron acceleration [Chaston et al, 2002] and 
in modifying F-region densities [Boehm et ai, 1990; Streltsov 
and Lotko, 2008]; however, their origin remains a subject of 
discussion. 

A particularly curious characteristic of such waves is their 
period, which was estimated for one event as 20 to 40 sec- 
onds [Karlsson et ai, 2004]. It is difficult to associate such 
periods with ionospheric timescales, such as trapping in the 
ionospheric Alfven resonator (lAR), because these tend to 
be about one second or less. Furthermore, waves are not 
confined to the lAR, being observed at altitudes up to 4 
Re. Attempts to associate waves with magnetospheric trap- 
ping likewise encounter difficulties, because typical magneto- 
spheric travel times are hundreds of seconds (e.g. ultra-low 
frequency waves in the Pc5 range) and waves occur near 
the boundary between open and closed field-lines. Thus 
attempts to give these waves a purely magnetospheric or 
purely ionospheric explanation are unsuccessful. 

One line of inquiry that has been promising is the premise 
that small-scale waves might be produced by self-consistent 
magnetosphere-ionosphere (M-I) coupling. This was inves- 
tigated by Streltsov and Lotko [2004]; Streltsov and Karls- 
son [2008] , whose numerical experiments demonstrated that 
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waves resembling observations are produced by a simu- 
lated system's response to large-scale field-aligned currents 
(FACs) via the ionospheric feedback mechanism (IFM). 

In their original paper, Streltsov and Lotko [2004] identi- 
fied instability as the most likely explanation for production 
of small-scale waves. This explanation is appealing because 
waves first appeared inside a region where the underlying E- 
region was depleted (low Pedersen conductance) and electric 
field was strong, and these conditions are known to favor 
growth of ionospheric feedback instability (IFI) [Atkinson, 
1970; Holzer and Sato, 1973; Sato, 1978; Lysak, 1991; Lysak 
and Song, 2002]. It was therefore proposed that the experi- 
ment might develop small-scale waves through amplification 
of initial seed perturbations by IFI. 

Some features of the waves observed in numerical experi- 
ments and the real magnetosphere do not fit well with an IFI 
explanation. Feedback instability requires a trapping region 
so that waves may be overrefiected from an underlying E- 
region many times. Trapping should therefore be expected 
to leave a signature on waves that have grown through IFI, 
e.g. by imprinting a period or magnetic-field aligned length 
scale on the resulting waves. The waves produced in the ex- 
periments by Streltsov and Lotko [2004] did not show such 
periodic structure along the field-aligned direction, and ob- 
served wave periods are an order of magnitude longer than 
typical lAR periods. 

The evolution of small-scale waves in the magnetosphere 
and ionosphere and their impact on auroral processes has 
also been studied by, e.g., Seyler [1990]; Chaston et al. 
[2002]; Genot et al. [2004]; Lysak and Song [2008] and 
references therein, often using highly sophisticated mod- 
els. These studies, however, all imposed either some short 
timescale or wavelength on the system through their driver 
or initial conditions. Therefore, although they offer valu- 
able insights into the evolution and impacts of small-scale 
Alfven waves, they did not address how the waves are ini- 
tially produced and how the system imposes timescales and 
wavelengths, which is the focus of this work. 

This paper aims to establish an alternative explanation 
for the origin of small-scale Alfven waves produced by non- 
linear M-I coupling, and to show how wave properties are 
determined. First, we demonstrate that feedback instabil- 
ity in not necessary to produce small-scale waves in down- 
ward current regions, by presenting a simplified numeri- 
cal experiment which excludes feedback instability while 
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retaining the feedback mechanism. Wc then ofTer a new 
interpretation in terms of ionospheric depletion, wavcbreak- 
ing (nonUnear steepening), magnetosphere-ionosphcrc waves 
and frequency-based phase mixing. This novel explana- 
tion draws on recent theoretical advances in the topic of 
magnetosphere-ionosphere coupling, which are linked for the 
first time. Using these, we are able to clarify the process 
by which M-I coupling can produce small-scale waves, and 
highlight formulas for wave periods and wavelengths. 

2. Modeling 

The first goal of this paper is to test the hypothesis that 
small-scale Alfvon waves seen in rmmerical M-I coupling ex- 
periments are produced by the action of ionospheric feed- 
back instability. 

IFI has two essential ingredients. First of all, the iono- 
spheric feedback mechanism is required to act in the pres- 
ence of an ambient electric field, which may be due to the 
presence of a large-scale current system. Under these con- 
ditions, small-scale Alfvcn waves, incident on the E-rcgion 
from above, can be ovcrrcflected if they have suitable prop- 
erties [Trakhtengerts and Feldstein, 1984; Lysak and Song, 
2002; Russell and Wright, 2012]. In these cases, the re- 
flected small-scale wave has greater amplitude than the in- 
cident wave, increased energy flux being accounted for by 
a reduction in the ionospheric heating caused by the ambi- 
ent electric field. Instability can follow if a second ingredi- 
ent is present: reflection of upgoing Alfven waves at some 
location above the E-region, e.g. from the conjugate iono- 
sphere or from gradients in Alfvcn speed below the peak 
at about 1 Rb in which case waves are partially trapped 
inside the lAR. If this second reflection produces a down- 
going wave whose amplitude is greater than the initial in- 
cident wave, then successive cycles of overreflection at the 
E-region and partial reflection at an overlying altitude cause 
waves to grow exponentially and a system seeded with low- 
amplitude perturbations becomes dominated by the fastest 
growing mode [Atkinson, 1970; Holzer and Sato, 1973; Sato, 
1978; Trakhtengerts and Feldstein, 1984; Lysak, 1991; Lysak 
and Song, 2002]. 

Since both IFM and wave trapping are required for IFI, 
the hypothesis that IFI is responsible for producing the 
small-scale waves seen by Streltsov and Lotko [2004] can be 
tested using a simplifled model that includes the essential 
elements of IFM, but which is designed to avoid reflections 
above the El-region, thus preventing IFI. 

To this end, we have investigated what happens when a 
sheet ionosphere is coupled to a uniform overlying plasma 
and the system driven with large-scale currents in the form 
of a large-scale incident Alfven wave. The sheet ionosphere 
description is mathematically valid provided the electric 
fleld skin depth in the ionosphere is larger than or com- 
parable to the ionospheric thickness, generally true for fre- 
quencies below 100 Hz [Lysak, 1991]. Matters are simplifled 
by using a 2D model, sketched in Figure 1, in which z is 
the vertical coordinate, x is an invariant horizontal direc- 
tion, and y completes the Cartesian system. Equilibrium 
magnetic fleld is assumed vertical (appropriate for polar lat- 
itudes and low altitudes) and the direction of the horizontal 
electric field defines the j/-direction. We arc interested in 
timescales much longer than cyclotron or plasma periods, 
so the plasma is described as a single fluid. Wave ampli- 
tudes are assumed small enough for nonlinear effects to be 
unimportant in the magnetospheric plasma region, because 
these are unlikely to play a role in producing small scales via 
IFM; the majority of current closure by Pedersen currents 
is assumed to occur in the sheet ionosphere, so that Peder- 
sen conductivity can be neglected in the overlying plasma; 



and we neglect dissipation due to Landau damping, viscos- 
ity and resistivity, which may affect the long term evolution 
of small-scale Alfven waves but arc unlikely to be important 
for their creation. The overlying plasma region is therefore 
described using a single fluid model in which dynamics can 
be described in terms of inertial Alfven waves. 

IFM is included by solving the ionospheric response to 
field aligned current through the height-integrated continu- 
ity equation: 



(1) 



where N is the height-integrated ionospheric plasma density, 
jz is the field-aligned current at the top of the sheet iono- 
sphere, e is the fundamental charge, a is an effective recom- 
bination coefficient [Brekke, 1997], and h is the ionospheric 
thickness. The source term aN'^/h represents background 
ionization, the loss term —aN^/h represents recombination 
and jz/e accounts for addition (removal) of electrons by up- 
ward (downward) FACs. 

Ohmic heating can change recombination rates, and this 
can be important for F-region modeling, where the low den- 
sities allow for significant changes on timescales as short as 
a few minutes [Saint Maurice and Torr, 1978; Zettergren 
et al, 2010; Zettergren and Semeter, 2012]. In the E-region, 
however, the much greater thermal mass means that a may 
be treated as constant for the timescales we study (several to 
at most tens of minutes). Heating of the F-region may affect 
the Alfven speed there, and we discuss the consequences of 
this in our discussion, but our primary conclusions regard- 
ing creation of small-scale waves and their properties are 
unaffected. 

The feedback loop of IFM is completed by imposing a 
boundary condition at the base of the overlying plasma that 
depends on ionospheric number density. This is derived us- 
ing current closure and takes the form. 



= {noeMp)NEy 



(2) 



for our 2D model, where Mp is the Pedersen mobility defined 
by Ep = eMpN with Ep the effective height-integrated 

Pedersen conductance. Mp is taken constant in our study. 

Magnetospheric equations (describing Alfven waves in a 
fully ionized plasma with electron inertia, under the assump- 
tions already stated) are as follows: 
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A derivation of these equations is given by Russell [2010] and 
they can also be found as a limiting case of the equations 
used by Lysak and Song [2008]. 

The whole system is solved numerically using an explicit 
leapfrog trapezoidal scheme and centered finite differences 
[Russell, 2010]. In equation (5), the parameter c is as- 
signed an artificial value that ensures the electron plasma 
frequency, LOp^ = c/Ag, is much greater than the angular fre- 
quency of the oscillations of interest (ensuring proper inertial 
Alfven wave behavior) but can also be resolved with a rea- 
sonable timestep (ensuring numerical stability). This tech- 
nique has been previously used by Lysak and Song [2008] 
and is discussed by Russell [2010]. Reflections from the 
upper boundary (which could potentially lead to a numeri- 
cal form of IFI) are avoided by positioning it at a distance 
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Figure 1. Geometry of simplified electrodynamic 
magnetospliere-ionospliere coupling model. 



that ensures waves cannot leave the subsection of the sim- 
ulation domain shown in our figures, reflect from the upper 
boundary and return to the subdomain within the simula- 
tion runtime. Thus, only part of the magnetospheric domain 
is shown in our figures. 

The resulting model is much simpler than others that 
have been used to describe the ionosphere-thermosphere sys- 
tem, e.g. Seyler [1990]; Dreher [1997]; Zhu et al. [2001]; 
Streltsov and Lotko [2004]; Sydorenko et al. [2008]; Chas- 
ton et al, [2011]; Zettergren and Semeter [2012], however, 
for the purposes of this paper it offers some important ad- 
vantages. First, it ensures that any small-scale waves pro- 
duced are created by IFM but not IFI, providing a rigorous 
test of whether or not IFI is required, and second, it allows 
the production of small-scale waves to be studied in isola- 
tion from their later development, which might otherwise 
obscure the details of their origin. Simple models have a 
strong track record for showcasing the essential aspects of 
a problem, and, if designed appropriately, they are valuable 
for developing new interpretations and identifying formulas 
that remain a useful guide when additional physics is later 
added. The task of adding greater complexity, and seeing 
how these fundamental behaviors are altered, is left for the 
future. 

It is at times useful to refer to simulations that do not in- 
clude electron inertia in the magnetospheric domain, partic- 
ularly to clarify how small scales are produced in the absence 
of IFI. These are performed using an approach proposed by 
Cran-McGreehin et al. [2007] which reduces the model to a 
single governing equation that is solved numerically using a 
leapfrog trapezoidal scheme with one-sided finite differences 
[Russell et al., 2010]. Apart from the neglect of electron 
inertia, the assumptions for this ideal MHD model are the 
same as for the electron-inertial single-ffuid model. 

3. Waves Produced in Absence of Instability 

To determine whether small-scale waves are produced by 
IFM in the absence of IFI, the electron-inertial M-I model 
described in Section 2 was driven with a system of large-scale 
currents and allowed to evolve. The currents were created 
by specifying a large-scale incident Alfven wave that shears 
the background magnetic field to produce one channel of up- 
ward FAC and one channel of downward FAC. (This form 
of driving is equivalent to that used by Streltsov and Lotko 
[2004].) The incident wave does not contain small scales and 
remains constant after a short ramping transient. Current 
systems like this may be produced by processes in the mag- 
netosphere, far from the ionosphere, but their origin is not 
considered here. Suffice to say that currents do form, and 
we examine the M-I system response. 

The evolution of our simulation is shown by Figure 2, 
which presents three snapshots of field-aligned current den- 
sity (normalized to the maximum current density of the in- 
cident driving wave, ji) and height-integrated ionospheric 



number density (normalized to the equilibrium value in the 
absence of field-aligned currents, A^e). The sheet ionosphere 
is positioned at 2 = 0, which corresponds to an altitude 
of approximately 110 km. Times are normalized by a de- 
pletion timescale, r = eN^/ji, this choice lending itself to 
the present study better than the advection timescale used 
for normalization by Russell and Wright [2012]. Horizontal 
distances are normalized by the width of the simulation do- 
main, j/o, and vertical distances are normalized by zq = ^at 
where va is the Alfven speed. Typically, zo 3> yo. 

Examining Figure 2, at t/r = 1.4, the large-scale current 
system is well established and ionospheric number density is 
evolving in response to the currents. One can see that A'^ has 
been increased in the upward current channel (centered on 
y/yo~ 0.3), as electrons are deposited there, and decreased 
in the downward current channel (centered on y/yo — 0.7), 
as electrons are removed there. 




0.4 0.6 

y/Yo 

t/r = 3.0 




0.4 0.6 

y/Yo 

Figure 2. Response of a uniform plasma with finite 
electron inertia and a sheet ionosphere to a large scale 
current. Color shows field-aligned current, j^, and line 
plots show height-integrated ionospheric plasma density, 
N. 
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At tjr = 2.2 the upward current channel has very nearly 
reached a new steady state (the maximum value of A*' is 
within 0.3% of its steady-state and converging slowly to it), 
however, the downward current channel and associated iono- 
spheric depletion region are in the process of broadening on 
the left-hand side, where they expand in the direction of the 
horizontal electric field. As they widen, becomes steep 
and a ripple begins to form around 2//1/0 = 0.55, at this 
time looking like a small tooth at the left-hand edge of the 
depleted region. 

At later times, e.g. tjr = 3.0, small-scale waves are 
clearly present. Small-scale changes in N, apparent between 
y/yo — 0.48 and y/yo = 0.58, correspond to small-scale 
changes in ionospheric reflectivity, so when the large-scale 
incident Alfvcn wave driving the system reflects from the 
ionosphere, small-scale upgoing Alfvcn waves arc formed 
(visible in the color plot of jz at t/r = 3.0 as alternating 
bands of blue and purple in the downward current channel 
with short horizontal wavelength). Those are incrtial Alfvcn 
waves and they carry small scales that originate at 2; = 
into the magnetosphere. Their phase-speed, transverse to 
the magnetic field, is in the direction of the large-scale elec- 
tric field (to the left) and their group velocity, transverse to 
the magnetic field, is oppositely directed (to the right), caus- 
ing them to spread out over the downward current channel 
with increasing altitude. 

These features axe similar to those seen at eaxly times in 
the experiments of Streltsov and Lotko [2004]. In particu- 
lar, boundary waves form just inside the downward current 
channel, at the edge adjacent to the upward current channel, 
and the same tooth-like feature is seen in both experiments 
(see their Figures 3 and 4). We therefore conclude that 
small-scale waves arc formed by the same process in both 
experiments, even though the experiments of Streltsov and 
Lotko [2004] included many more aspects of the M-I sys- 
tem. This reassures us that although the present model is 
simpler, it does include the essential elements we wish to 
understand. 

These results demonstrate that small-scale waves can be 
produced by nonlinear M-I coupling in the absence of insta- 
bility, claxifying the mechanism of their origin. In particular, 
they do not arise through amplification of noise by IFI as 
suggested previously. We now outline a new interpretation 
that proposes how small-scale waves may be formed and 
what their properties will be. 

4. Small Scales from Ionospheric Depletion 



Formation of small scales appears to involve ionospheric 
depletion and horizontal expansion of the depleted region 
and associated downward current channel. It is therefore 
worthwhile to examine the conditions under which broaden- 
ing occurs and the manner in which it happens, with par- 
ticular focus on the generation of small scales. 

The theory of ionospheric depletion and cavity /current 
broadening has previously been studied by Doe et al. [1995 
BUxt and Brekke [1996]; Karlsson and Marklund [1998 
Cran-McGreehin et al. [2007]; Russell et al. [2010]; Zetter- 
gren and Semeter [2012] and we quote their work for the 
condition for broadening. These authors have shown that 
the coupled M-I system has two distinct typos of evolution, 
depending on the maxirrmm strength of downward current 
density. Examining equation (1), background ionization can 
balance downward FAC provided that \jz\ < jc, where 



h 



(7) 



For an ionosphere that is initially highly conducting, the 
maximum current density at its top is approximately 2ji, 
where ji = msLx(ahs{d(SBi) / dy)) / no is the maximum paral- 
lel current density of the incident Alfven wave. If \2ji\ < jc, 
then the coupled system evolves smoothly to a new partially- 
depleted steady state [Cran-McGreehin et al, 2007; Russell 
et at, 2010]. 

If > jc, then dynamics arc more dramatic. Fig- 

ure 3 shows a simulated system response to a downward 
current density that exceeds the critical value and is there- 
fore unsustainable. Starting with an initially uniform iono- 
sphere, downward FAC strongly depletes ionospheric plasma 
to such an extent that ionospheric reflectivity is significantly 
reduced, decreasing downward current density; hence, the 
weakened downward FAC is forced to broaden to close the 
upward FAC. Here, electron inertia is neglected and as the 
downward FAC broadens, a moving density discontinuity 
forms in the ionosphere, corresponding to a moving current- 
sheet in the magnetosphere. 

The ideal simulation presented in Figure 3 uses an ini- 
tial condition with Ep/Sa = 100, for which the features of 
broadening are very clear. Our inertial model is not capa- 
ble of solving under these conditions because deep depletion 
associated with high values of Ep/E^ can trigger numerical 
instability through our boundary condition. Consequently, 
the sirrmlation shown in Figure 2 uses Ep/Ea = 10 for sta- 
bility, while the simulation shown in Figure 2 uses a higher 
value for improved clarity. This means that these figures 
should not be compared quantitatively, although a qualita- 
tive comparison is very valuable. 

The aspects of broadening outlined so far have been de- 
scribed previously by Cran-McGreehin et al. [2007], how- 
ever recent theoretical developments allow us to add some 
new comments about the broadening process. First, it is 
interesting to consider broadening in the context of recent 
knowledge that the governing equations for IFM reduce to 
an advection equation in the absence of electron inertia [Rus- 
sell and Wright, 2012]. Thus, dynamics arc governed by a 
characteristic advection speed, which, for a ID sheet iono- 
sphere, produces motion in the direction of the electric field 
at a speed 



vmi 



MpEy 

1 + Ep/E^ 



(8) 



Consequently, the ionospheric response to an incident Alfven 
wave depends on the strength of the current density the 
large-scale incident wave creates in relation to this thresh- 
old. 



[Russell and Wright, 2012]. Preexisting arguments that ex- 
plain broadening only by need to provide current closure 
(such as that outlined above) do not lead one to expect 
asymmetry. However, since the advective properties only 
permit motion in the direction of the electric field, it follows 
that broadening can only happen at an edge where the elec- 
tric field points away from the depleted region. This explains 
why the depleted region and downward current channel may 
broaden only on one side. 

The discontinuity itself moves at the M-I advection speed 
derived for discontinuities by Russell and Wright [2012], 
which is the geometric mean of the values of vmi to either 
side of the discontinuity. This gives the rate of expansion 
and is applied here to ionospheric depletion for the first time. 

It is also informative to compare a single snapshot from 
Figure 3 to the self-consistent M-I steady states obtained 
by Russell et al. [2010]. Such a comparison in provided in 
Figure 4. From this, one sees that the discontinuity steps 
between an "upper" steady state (computed under the as- 
sumption Ep S> Ea and plotted in blue) and a "lower" 
steady state (computed under the assumption A'^^ <C 
and plotted in green). When the discontinuity first forms, it 
does so where the upper and lower steady states cross inside 
the downward current channel (at y/yo = 0.59). Thereafter, 



RUSSELL ET AL.: PRODUCTION OF SMALL-SCALE ALFVEN WAVES 



X - 5 



the discontinuity sweeps to the left across the region where 
both steady states are valid and stops where they cross again 
(at y/yo = 0.49). Thus, these solutions indicate the early 
and final widths of the depleted region, as well as the height 
of the discontinuity at any position. 

The combination of advective properties with the exis- 
tence of two steady state solutions, leads us to propose 
that broadening is essentially an M-I wavebreaking (non- 
linear steepening) process. A useful visualization of ad- 
vection, sketched in Figure 5, is to imagine the motion of 
points over time in (y, N) space under the action of ad- 
vection in y {dy/dt = vmi) and changes in due to the 
convective derivative, dN/dt. In an infinitesimal time St, a 




Figure 3. Evolution with incident downward current 
satisfying \2ji\ > jc and no electron inertia, showing for- 
mation of a discontinuity in ionospheric number density 
and a corresponding current sheet. Snapshots show iono- 
spheric plasma density (bottom) and field-aligned current 
at the top of the ionosphere (top). 




y/Yo 



Figure 4. Snapshot of density discontinuity formed 
in absence of electron inertia (red), upper steady state 
(blue) and lower steady state (green). 



point initially at (y,N) moves to (y -I- Sy,N + SN), with 
Sy = VMi{y,N)5t and 5N = {dN/dt)5t. Thus, in dynamic 
evolution (panel (a) of Figure 5) , N{y,t) evolves as the curve 
is carried by the motion of the points it threads. In this con- 
text, a steady state solution (panel (b) of Figure 5) is a tra- 
jectory common to all points lying on it, hence, although 
points move along the steady state, the curve threading 
them does not change, giving dN/dt = 0. 

Referring to Figure 3, our simulation commences with dy- 
namic evolution, as points starting on the initial condition 
move towards the nearest steady state. In the upward cur- 
rent channel and the at edges of downward current channel, 
this is the upper steady state. In the center of the downward 
current channel, it is the lower steady state. At the right 
hand edge of downward channel, the upper and lower steady 
states break down at a common location and the true steady 
state is a matching between these which can be obtained us- 
ing the asymptotic method of Russell et al. [2010]. At the 
left hand edge of downward channel, the upper steady state 
is encountered first, so points in {y, N) space pause here 
although there is also a valid lower steady state at lower A'^. 

Once points have moved on to a steady state curve, they 
move along it, with motion in y determined by the advec- 
tion speed. The vital property responsible for wavebreaking 
is this: examining equation 8, a point at {y, N) moves more 
rapidly than a point at the same y but greater A'^ (for which 
Ep/Ea is higher). Thus, in the region where there are two 
valid steady states, points following the upper steady state 
move more slowly than, and are overtaken by, points moving 
along the lower steady state (panel (c) of Figure 5). When 
points catch up in y, a discontinuity forms to preserve a sin- 
gle valued physical solution. This wavebreaking (nonlinear 
steepening) picture presents a new conceptual view of the 
broadening that occurs in response to ionospheric depletion, 
connecting broadening with the advective properties of M-I 
coupling and the existence of two steady state solutions. 

We conclude from this section that the coupled M-I sys- 
tem produces small horizontal scales in response to strong 
downward FAC that exceeds the critical threshold jc, and 
that it does this through ionospheric depletion and broad- 
ening of downward FAC, which can be considered as a wave- 
breaking (nonlinear steepening) process, that occurs because 
of the existence of two steady state solutions. 



, (a) Dynamic evolution . , (b) Steady state 
L N 
A A 




(c) Nonlinear steepening 
'a from two steady states 




Figure 5. Trajectories of points in {y, N) space under 
advection. (a) Motion of points during dynamic evolu- 
tion, giving dN/dt ^ 0. (b) Motion of points along a 
steady state curve, maintaining dN/dt — 0. (c) Where 
two steady states exist, points move more rapidly along 
the lower steady state, so a moving discontinuity forms 
to preserve a single valued solution. 
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5. Wave Properties and Evolution 

In reality, gradient length scales do not collapse to zero 
because small-scale physics prevents this. For the present 
M-I coupling work, the appropriate modification is to in- 
clude electron inertia in the magnetosphere. 

The changes due to electron inertia can be seen by com- 
paring the qualitative features of Figure 2 (finite electron in- 
ertial length, \e/yo ~ 0.025) with those of Figure 3 (Ae — 0). 
At early times, length scales seen in Figure 2 arc much 
greater than Ae, so the simulation follows the ideal MHD 
evolution described in Section 4. When wavebreaking oc- 
curs, a steep transition appears and moves to the left, but it 
has finite width comparable to Ae instead of being a discon- 
tinuity. This broadening front is trailed by an undershoot, 
which later develops into a train of undershoots and over- 
shoots. Time evolution at a fixed location is gradual until 
the broadening front passes, its movement causing a rapid 
decrease in TV that produces oscillations. Such effects are 
to be expected by analogy to MHD shocks, which exhibit 
very similar features when the dominant small-scale physics 
is dispersive. 

The wake of undershoots and overshoots that trails the 
broadening front in N is of great interest because these dis- 
turbances drive the upgoing inertial Alfven waves whose ori- 
gin and properties we aim to explain When the incident 

large-scale wave driving the system reflects from an iono- 
sphere with small-scale variations in ionospheric reflectivity 
(due to small scale variations in A^), these variations impose 
their scale and period on the reflected wave, which carries 
these properties into the magnetosphere. 

Wave-like behaviors resulting from self-consistent elec- 
trodynamic M-I coupling were recently described by 
Russell and Wright [2012], who suggested the name 
"magnetosphere-ionosphere (M-I) waves" for the phenom- 
ena and also identified nonlinear steepening as a means of ac- 
cessing small scales. The properties of M-I waves are shared 
by the small-scale disturbances produced on the boundary 
in the present study, and it is useful to refer to them for in- 
terpretation. The current work docs, however, advance our 
understanding of these waves in three important ways: (i) it 
demonstrates that small-scale M-I waves can be be produced 
from a simple large-scale driver by the depletion mechanism 
outlined above; (ii) the spatially varying steady state intro- 
duces M-I wave phase mixing (discussed below) that reduces 
length scales further and was not considered by Russell and 
Wright [2012]; and (iii) we introduce a simple formula for 
the largest inertial wavelength produced by wavebreaking. 
We also note that Russell and Wright [2012] did not apply 
their results to the generation of small-scale Alfven waves 
launched into the magnetosphere following ionospheric de- 
pletion, which is the focus of this paper. 

According to the theory of M-I waves, as electron inertial 
becomes significant, the advection that occurs for Ae = 
is altered, giving way to oscillation at a characteristic fre- 
quency for wavelengths comparable to or less than Ae (see 
Figure 4 of Russell and Wright [2012]). This frequency is 
given by 



fMI = 



MpEy 



27rAe (Ep/E/ 



(9) 



Meanwhile, wave amplitudes decay on the ionospheric re- 
combination time, 



'^decay — 



2aN' 



(10) 



The time interval is chosen to focus on the evolution of the 
oscillations and not the gradual evolution and rapid fall in 
N that precedes them. Vertical dashed lines are separated 
by tmi/t where tmi = ^1 fhii is the period expected for 
M-I waves, and the amplitude envelope corresponds to ex- 
ponential decay with an e- folding time of Tdecay Both /m/ 
and Tdecay are evaluated using steady state values of Ey 
and N about which the waves oscillate. The formulas for 
/m/ and Tdecay (equations (9) and (10) respectively) are an 
excellent match to the oscillations, even though background 
quantities vary in y and the formulas were derived using nor- 
mal mode analysis, effectively looking for local linear modes 
about a global equilibrium. 

The high level of agreement gives us confldence in con- 
necting the properties of inertial M-I waves with small-scale 
waves produced in response to strong downward current. 
This link, which is new, is significant because: (i) it provides 
formulas for the properties of small-scale waves produced by 
the mechanism described in Section 4; and (ii) it identifies 
a natural process capable of creating M-I waves, which fur- 
thers their study. Equation (9) is easily used and gives fre- 
quencies that are likely to dominate a low quality lAR, (i.e. 
an lAR for which Alfven speed gradients do not produce 
substantial rcficctivity for the wave frequencies under con- 
sideration, implying long growth times for IFI and allowing 
substantial transmission of upgoing waves to the magneto- 
sphere). A high quality lAR with large IFI growth rate 
could also exhibit the conventional fastest growing mode 
frequency, computation of which requires more cumbersome 
numerical simulations [Streltsov et ai, 2010]. 

It remains to determine the wavelength of the M-I waves 
on the boundary. Figure 7 shows a close-up of the pertur- 
bation in N at two different times, with the electron iner- 
tial length indicated. The waves clearly have a wavelength 
(measured along the boundary) comparable to Ae and their 
group velocity is negligible. Importantly, wavelengths be- 
come shorter over time. 

The change in wavelength of M-I waves is a novel mani- 
festation of frequency-based phase mixing. Since the back- 
ground values of Ey and N vary with y, so too does /mi 
(see equation 9), as is apparent from comparing oscillations 
at y/yo ~ 0.52 and y/yo = 0.54 (Figure 6). It follows 
that M-I waves, formed on the boundary with wavelengths 
slightly larger than the electron inertial length at the base of 
the magnetosphere, phase mix to ever smaller length scales. 
Frequency-based phase mixing has been well studied for field 
line resonances [e.g. Mann et al., 1995; Rankin et al., 2005; 
Russell and Wright. 2010] but we apply it for the first time 
to M-I waves, drawing on the result. 



{t\dfMildy\y 



(11) 



We find the wavelength prior to phase mixing (the length 
scale immediately trailing the broadening front) can be es- 
timated by equating the speed of an equivalent discontinu- 
ity with the phase speed of the trailing disturbances. This 
equality ensures that the broadening front and the under- 
shoot immediately behind it move together as a single struc- 
ture. 

In the ideal limit, a discontinuity moves at a speed given 

by 



Vd = VMiia)vMi{b), 



(12) 



Figure 6 plots the relative perturbations of N at y/yo = 
0.52 and y/yo = 0.54 between times t/r = 2.8 and t/r = 5.8. 



where "a" ("b") indicates a function or value is evaluated 
ahead of (behind) the discontinuity and vmi is given by 
equation (8) [Russell and Wright, 2012]. When electron in- 
ertia is included, the broadening front is smoothed from a 
discontinuity, however equation (12) may be used as a first 
approximation. Ripples behind the broadening front evolve 
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as M-I waves, which have a phase speed given by Russell 
and Wright [2012] as 



MpEy 



l + (Ep/SA)Vl + fciAi 



(13) 



Equating Vp (equation (12)) and Vd (equation (13)), and 
making use of the relationship Ey — (1 + r)Ei , where 
r = (1 — Ep/Ea)/(1 + Ep/Ea) is the electric field reflection 
coefficient and the incident electric field Ei is approximately 
the same ahead and behind the broadening front, the wave- 
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Figure 6. Perturbation in height-integrated ionospheric 
number density relative to steady state value, at y/yo ~ 
0.54 (top) and y/yo = 0.52 (bottom). For each panel, 
vertical dashed lines are separated by tmi/t, and the 
dashed envelope is calculated for exponential decay at the 
recombinative e-folding time Tdecay/T, where these times 
are computed from steady state values at each location. 



0.10 



0.05 



CO 
CO 



> 0.00 



-0.05 



-0.10 





: th = 3.3 : 




'\ : t/T = 4.4 : 







length of disturbances produced by broadening is 



Xy — 



1 



(14) 



being determined solely by Ae and the ratio of number den- 
sities (equivalently Pedersen conductances) ahead of and be- 
hind the broadening front. 

Equation (14) can be tested against the inertial sim- 
ulation presented in Figure 2. At time t/r — 2.2, the 
wavelength behind the broadening front is estimated from 
perturbed number density as twice the distance between 
the local minimum immediately following the front and 
the following local maximum. This gives Ay/Ae = 2.0. 
The number density behind the front is easily identified as 
N{b)/Ne ~ 0.312 (the local minimum). For this simulation, 
the value ahead of the front is less precisely determined be- 
cause there is some ambiguity about where the transition 
begins. We therefore give the range 0.8 < Na/Ne < 1.1, for 
which equation (14) returns 1.86 < Ay/Ae < 2.66, in good 
agreement with Ay/Ae = 2.0. 

This wavelength analysis can also be applied when an 
initially large-scale M-I wave steepens nonlinearly to form 
ripples behind a steep gradient, as seen in Figure 2 of Lysak 
and Song [2002] and Figure 6 of Russell and Wright [2012]. 
In the latter case, small scales have formed hy t/r = 4, at 
which time the density ratio across the acting discontinuity 
is N{a)/N{b) = 1.284 (the value for N{a) is precisely deter- 
mined from a local maximum ahead of the main transition). 
Therefore, equation (14) returns Xy/Xe = 7.8, which agrees 
exactly with the wavelength measured from crossings of A'^ 
with Ne- This agreement gives further confidence in equa- 
tion (14) and demonstrates its wider applicability. Equation 
(14) also explains why Xy/Xe is nearly four times larger for 
the wavebreaking simulation of Russell and Wright [2012] 
than for the ionospheric depletion scenario studied in the 
present work: depletion leads to much greater N{a)/N{b), 
so Xy/Xe is scaled accordingly. 

Thus, production of small scales is a two stage process. 
First, strong downward FAC causes a rapid collapse of hori- 
zontal scales, populating the M-I boundary with waves that 
have a horizontal scale of approximately the electron iner- 
tial length for the low-altitude magnetosphere, which drive 
upgoing inertial Alfven waves. Thereafter, frequency-based 
phase mixing of boundary waves continuously reduces hori- 
zontal scales further. 



6. Discussion 
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Figure 7. Perturbation in height-integrated ionospheric 
number density relative to steady state value, at t/r = 
3.3 (red) and t/r = 4.4 (blue) with dashed lines showing 
the electron inertial length. 



The present work shows that small-scale waves are pro- 
duced by the self-consistent M-I response to strong down- 
ward current, even when amplification due to IFI is ex- 
cluded. This suggests that instability is not responsible for 
producing these waves, although nonlinear M-I coupling is 
most definitely at work. We propose that small scales are ac- 
tively created by a nonlinear steepening/wavebreaking pro- 
cess that follows ionospheric depletion, which is also respon- 
sible for expanding the ionospheric cavity and associated 
downward current channel in the direction of the electric 
field. In this view, broadening rapidly collapses the short- 
est horizontal gradient length scale to approximately the 
O.BOelectron inertial length for the low-altitude magnetosphere, 
the exact scale determined by the ratio of ionospheric num- 
ber densities (equivalently Pedersen conductances) to either 
side of the broadening front. After inertial waves have been 
formed, length scales are further diminished by frequency- 
based phase mixing of M-I boundary waves. This interpreta- 
tion contrasts with a previously held view that ionospheric 
depletion simply creates favorable conditions for IFI, and 
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leads to different conclusions about resulting wave proper- 
ties. 

In the new interpretation, broadening produces boundary 
waves with the properties of inertial M-I waves, which oscil- 
late at a frequency determined by magnetospheric and iono- 
spheric properties. This frequency provides a central test as 
to whether this mechanism may explain previous observa- 
tions of small-scale waves in downward current channels. 

We shall estimate this frequency under conditions appro- 
priate to the observations of Karisson et al, [2004] , for which 
wave periods were reported as 20 to 40 seconds. At the time 
of observation, the Cluster satellites were passing over the 
nighttime winter southern hemisphere (magnetic footpoints 
at about 20:00 MLT and 70° S geomagnetic latitude at about 
05:30 UT on 19 May 2002). Under these conditions, E and 
F regions are normally well defined, so we will assume the 
sheet ionosphere in our model represents the E-region, and 
properties corresponding to the low altitude rnagnctosphcre 
will be evaluated at electron density minimum between the 
E and F regions. Using typical values, we take magnetic 
field strength B — 5.5 x 10^^ T, a low-altitude magneto- 



spheric electron density of n„ 



5 X 10 m and an av- 



erage ion mass mi = 30mp where nip is the proton mass (a 
good approximation for a plasma dominated by NO"*" fol- 
lowed by O^). This gives an electron inertial length of 240 
m and an Alfvcu conductance, Ea = 0.081 mho. Putting 
Mp = 10* m^s"^V"\ Sp = 0.6 mho and Ey = 0.1 V (for 
quantities in the depleted E-region), equation (9) returns 
a period of 11 seconds. This is very close to the observed 
period of 20 to 40 seconds and certainly within the uncer- 
tainty resulting from plausible deviations from our choice of 
typical values (which in any case will vary across the region 
populated with M-I waves). 

We also note that the F-region commonly becomes de- 
pleted in association with El-region depletion, through mo- 
tion of plasma along the magnetic field and enhanced re- 
combination due to Ohmic heating (as discussed in Section 
2). Since /mi oc Umag this effect can easily increase periods 
by a factor 2 or 3 to agree exactly with those reported by 
Karisson et al. [2004]. In fact, in the real system, we spec- 
ulate that one may find wave periods becoming longer over 
several minutes as the F-region depletes. 

Prom the frequency of M-I waves and the perpendicular 
wavelength at the top of the ionosphere, one can deduce the 
properties of upgoing inertial Alfven waves produced by the 
ionospheric depletion mechanism. This provides useful tests 
for identifying observations with this theory. We start by 
noting that since M-I waves have a phase speed in the di- 
rection of the ionospheric electric field (equation (13)), the 
phase-speed of the magnetospheric waves should also have a 
component parallel to the electric field (never anti-parallel). 

Wavelengths in the rnagnctosphcre arc now considered. 
As an inertial Alfven wave propagates out into the magne- 
tosphere, its wavelength perpendicular to the magnetic field 
may increase with height due to expansion of the flux tube, 
and also due to spreading caused by the perpendicular group 
velocity of the waves (seen in Figure 2). Spreading due to 
group velocity can be estimated from the ray angles at the 
edges of the upgoing wavepacket. Using the dispersion rela- 
tion for inertial Alfven waves, 



knVA 



(15) 



and the group speed expressions Vg^± = doj/dk±, Vg^\\ = 
duj/dk\\, one finds 



Vg,± _ 



k±\e 



(16) 



Setting uj = 2tt fMi under a WKB assumption and putting 
tyi)ical values for Ae and va, equation (16) returns values for 
iig^x/fgjl on the order of 10"'* or less (the term in brackets 
is a factor between zero and one). Thus, rays are close to 
field aligned, so that expansion of the magnetic field with 
height is likely to dominate over spreading due to transverse 
group velocity. 

If transverse Alfven speed gradients are present in the 
rnagnctosphcre then increases to transverse wavelengths of 
upgoing inertial Alfven waves may be partially countered by 
velocity-based phase mixing, which acts to reduce transverse 
wavelengths [Genot et al, 2004; Lysak and Song, 2008]. 
These arguments suggest that small-scale waves observed 
in the magnetosphere should have transverse wavelengths 
that map to an ionospheric wavelength similar to or less 
than the electron inertial length in the low altitude magne- 
tosphere if they are produced by the depletion mechanism 
we have identified. This should be verified in future by more 
detailed modeling designed to capture changes to the upgo- 
ing waves with height, at which stage damping processes, 
such as Landau damping, may also be considered. 

The models used in this paper have treated the M-I sys- 
tem as having a uniform magnetosphere overlying a sheet 
ionosphere. Upgoing Alfven waves were therefore able to 
travel unhindered to the outer magnetosphere. In reality, 
the F-region produces a steep gradient of Alfven speed, 
which may partially refiect Alfven waves. Thus, Alfven 
waves may be partially trapped between the E and F re- 
gions, undergoing multiple rcfiections in a volume known as 
the ionospheric Alfven resonator (lAR) [Lysak, 1991]. 

To sec how F-rcgion trapping affects production and evo- 
lution of small-scale Alfven waves, it is useful to compare our 
findings with the simulations of Streltsov and Lotko [2004], 
e.g. their Figures 3 and 4. Their study included a non- 
uniform magnetosphere, F-region, 2D dipole field and 2-fiuid 
effects, and its output may be compared to the simulations 
presented earlier in Sections 3 and 4. 

Many features of the evolution are shared with our mod- 
els: FACs modify E-region plasma density, downward FAC is 
unsustainable, depleting the E-region and leading to broad- 
ening, and this, combined with electron inertia, populates a 
part of the M-I boundary with M-I waves that drive upgo- 
ing inertial Alfven waves that reach into the magnetosphere. 
These similarities suggest that our new interpretation, al- 
though founded on comparatively simple M-I models, does 
indeed hold when additional physics is included. 

The primary difference is that the simulations of Streltsov 
and Lotko [2004] show small-scale ionospheric disturbances 
gradually spreading over time, coming to occupy a substan- 
tial part of the downward current channel. This can be 
understood as follows: waves produced by M-I interactions 
have a component of phase-speed in the direction of the 
background electric field; therefore, inertial Alfven waves 
that are driven by M-I interactions have a component of 
group velocity in the direction opposite to the horizontal 
electric field. Consequently, during multiple reflections in- 
side the lAR, inertial Alfven waves will slowly spread across 
the downward FAC, modifying the ionosphere and extending 
the part of the boundary populated with M-I waves. 

We also note that IFI is able to act when F-rcgion reflec- 
tions are present. Thus, although IFI may not be respon- 
sible for the initial creation of small-scale waves, it could 
potentially affect their later evolution. The effects of IFI 
will depend on the steepness of the F-region Alfven speed 
gradient producing the trapping, and also the recombinative 
decay time in the E-region. If decay due to E-region effects 
is more rapid than growth due to IFI, we speculate that IFI 
will simply prolong the lifetime of small-scale waves by coun- 
teracting their decay. If IFI is strong enough to overcome 
decay processes, then it may significantly amplify small-scale 
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waves that have been produced by the processes described 
in this paper. 

More sophisticated modeling of the phenomena described 
in this paper is intended in future work. For example, 
our use of a height-integrated ionosphere, sufficient to cap- 
ture the most fundamental aspects of the present problem, 
nonetheless obscures issues such as at what altitude the iono- 
spheric density depletions are most extreme and the height 
profile of the disturbances that form the ionospheric portion 
of the M-I waves. Progressing to height-resolved models 
that include these features is therefore an important goal 
for future theoretical work. 

Other obvious developments axe to remove the 2D as- 
sumption made in this paper and to allow nonlinear effects 
in the magnetosphere. We speculate that events in such sim- 
ulations (using an initially 2D setup as in the present paper) 
would follow a similar initial evolution to those shown here. 
Once small-scale waves form, however, these narrow sheets 
will be subject to instabilities, including shear and tear- 
ing mode instabilities [Seyler, 1990; Chaston et al, 2011]. 
These may break up the 2D structures into filaments with 
azimuthal structure, for which M-I wave properties are less 
clear. 

It is also desirable that a detailed comparison to obser- 
vations be made. An ideal study would combine measure- 
ments of small-scale waves inside a large-scale downward 
current channel with simultaneous data from the E and F 
regions at the magnetic footpoint. The tests that would 
associate these waves with the processes described in this 
paper are: close agreement between measured (non Doppler 
shifted) frequencies and /m/ computed from measurements; 
a transverse length scale that maps to approximately the 
electron inertial length at the base of the magnetosphere; 
Vphase ■ E± > 0; and an upward Poynting flux carried by the 
small-scale waves, although care is needed because the total 
Poynting flux (due to small-scale waves and the large-scale 
driver) is downwards. 

To conclude this paper, we comment that understanding 
small-scale waves produced by M-I interactions could be sig- 
nificant for a number of disciplines. One example is predic- 
tion of ionospheric densities, which impact on communicar 
tions, GPS and satellite drag. Changes in ionospheric num- 
ber density, caused by FAC, are obvious in our simulations, 
but there is an additional aspect that we have not consid- 
ered: if small-scale Alfven waves attain sufficient amplitude, 
then they will exert a ponderomotive force on plasma above 
the E-region, moving it along the background magnetic field 
[Boehm et al., 1990; Streltsov and Lotko, 2008]. Our results 
clarify a mechanism by which Alfven waves performing this 
role may be produced, and constrain their frequency and 
wavelengths. 

The processes described here may also explain features 
of auroral dynamics. Further work is needed to model the 
auroral signatures of these waves, but the evolving current 
patterns are a guide to what may be expected. If one asso- 
ciates strong upward currents with optical auroral emission, 
one can imagine a preexisting auroral arc brightening as 
the large-scale current supplying it with electrons intensifies. 
Were current density in the corresponding downward FAC 
to increase beyond the critical threshold, the dark region 
would expand into the preexisting arc, while new narrow au- 
roral arcs, associated with the processes we have described, 
formed at its edge. As we have seen, the narrow arcs would 
move in the direction of the large-scale electric field (with 
the phase motion of the waves) but could spread in the op- 
posite direction (due to reflections of inertial Alfven waves 
in the lAR). Such a scenario is reminiscent of observations 
by Serneter et al. [2008] and would be an interesting topic 
for future investigation. 
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